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Abstract 

We demonstrate a method for calculating the neutral i?-meson decay constants and mixing 
matrix elements in unquenched lattice QCD with domain-wall light quarks and static ^-quarks. 
Our computation is performed on the "2+1" flavor gauge configurations generated by the RBC 
and UKQCD Collaborations with a lattice spacing of a ~ 0.11 fm (a~^ = 1.729 GeV) and a 
lattice spatial volume of approximately (1.8 fm)^. We simulate at three different light sea quark 
masses with pion masses down to approximately 430 MeV, and extrapolate to the physical quark 
masses using a phenomenologically-motivated fit function based on next-to-leading order heavy- 
light meson SU (2) chiral perturbation theory. For the 6-quarks, we use an improved formulation 
of the Eichten-Hill action with static link-smearing to increase the signal-to-noise ratio. We also 
improve the heavy- light axial current used to compute the i?-meson decay constant to 0{aspa) 
using one-loop lattice perturbation theory. We present initial results for the SU (3)-breaking ratios 
fBs/fSi and ^ = fBsyBBs/fBdyBBd^ thereby demonstrating the viability of the method. For the 
ratio of decay constants, we find fss/ fsa — 1.15(12) and for the ratio of mixing matrix elements, 
we find ^ = 1.13(12), where in both cases the errors reflect the combined statistical and systematic 
uncertainties, including an estimate of the size of neglected O^l/mj,) effects. 

PACS numbers: 12.38.Gc, 12.15.Hh, 14.40.Nd 
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I. INTRODUCTION 



Neutral S-meson mixing is a sensitive probe of quark flavor- changing interactions. When 
combined with experimental measurements of the and oscillation frequencies, pre- 
cise QCD determinations of the B^ and i?^-mixing hadronic matrix elements allow clean 
determinations of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements \Vtd\ and 
\Vts\ with all sources of systematic uncertainty under control. In the Standard Model, the 
mass-difference of the neutral i?-meson mass eigenstates Arriq (often called the oscillation 
frequency) is given by [3|, |j] 

Am, = Y^^^ KVt,\'So{x,)vBM,, (1) 

where q G {d, s}. Both the Inami-Lim function So{xt) with Xt = mf/m'^ jsl and the QCD 
coefficient r]B can be computed in perturbation theory 
AB = 2 mixing matrix element. 



The hadronic contribution to the 



m, = (bI 



[67^(1 -75)g][HXi-75)g]|5;), (2) 



must be calculated nonperturbatively from first principles using lattice QCD. 
The hadronic matrix element A^, is conventionally parametrized as 

M, = lmlJlBs,, (3) 

where m^^ is the mass of the 5- meson, fsg is the 5- meson decay constant and Bb^ is 
the 5-meson bag parameter. Many statistical and systematic uncertainties from lattice 
calculations cancel in the S'f/(3)-breaking ratio Aig/Md, which would be one in the limit 
rrid — mg. It is therefore convenient and conventional to introduce the quantity 



jBdAZ-DBa 

Given improved lattice calculations of the nonperturbative factor C,, recent experimental 



measurements of the oscillation frequencies Arrid and Arris to about 1% accuracy 



m 



now 



allow the possibility of precisely determining the ratio of the CKM matrix elements | Vtd|/|V(<j| 
(see, for example, Ref. jlo||). 



The ratio |Vtd|/|VJ<(| constrains the apex of the CKM unitarity triangle ll|, ll2j. It is likely 



that new physics would give rise to new quark-flavor changing interactions and additional 



CP-violating phases; these would manifest themselves as apparent inconsistencies among 
different measurements of quantities which should be identical within the standard CKM 
picture. Thus a precise determination of the ratio ^ will help to constrain physics beyond 
the Standard Model. Furthermore, possible indications of new physics in B^-mixing at the 
~2.7-a level Q make the lattice calculation of B -meson mixing parameters timely. 

Recently, the HPQCD Collaboration published the first unquenched determination of ^ 
with an accuracy of 2.6% [l^, and the Fermilab Lattice and MILC Collaborations expect to 



have a result with similar errors soon [15|. The HPQCD calculation employs a nonrelativistic 
QCD (NRQCD) action for the heavy 6-quark [igI, while the Fermilab/MILC calculation uses 
the relativistic "Fermilab" action for the 6-quark [l7|. Both of these computations, however, 
rely on the same "2+1" flavor asqtad-improved staggered ensembles generated by the MILC 
Collaboration , which include the effects of two degenerate light quarks and one heavier 
close-to-strange quark in the sea sector. 

For such a phenomenologically important quantity as ^, it is valuable to have an inde- 
pendent crosscheck using different formulations of the lattice action for both the light and 
heavy quarks. Our calculation employs the 2+1 flavor dynamical domain- wall ensembles 
generated by the RBC and UKQCD Collaborations with a lattice spacing of a ~ 0.11 fm 



(a~^ = 1.729 GeV) [19|]. The use of domain-wall fermions 2^22| has the advantage over 
other light-quark formulations that the chiral perturbation theory expressions needed to ex- 
trapolate domain-wall lattice results to the physical u- and d-quark masses are closer to the 
continuum forms and have fewer parameters than in the Wilson or staggered cases 23|, |24 |. 
We compute the 6-quarks in the static limit (m^ — )■ oo), which leads to correlation func- 



tions that are noisier than those with propagating 6-quarks such as in the Fermi 
or NRQCD actio ns jl6| . We therefore use the static-quark formulation of Refs. 

or HYP 



with either APE 



27 



ab [l7| 

26| 



25 



smearing of the static-quark gauge links to increase 
the signal-to-noise ratio and reduce scaling violations (for some quantities) as compared to 
the Eichten-Hill action jsQl. Furthermore, the approximate chiral symmetry of the domain- 
wall action combined with the spin symmetry of the static action simplifies the lattice-to- 
continuum operator matching as compared to the Wilson case by reducing the number of 



additional lattice operators which appear 3l|, |32|. The results of this work extend an earlier 



study with two flavors of dynamical quarks and heavier light-quark masses in Ref. 



33|. 



The primary purpose of this paper is to demonstrate the viability of our method for 



4 



computing the S-meson decay constants and Ai? = 2 mixing matrix elements. We therefore 
use the small-volume (16^) ensembles with only one lattice spacing and have relatively 
heavy light-quark masses and limited statistics. A novel feature of this work is the use of 
SU{2) heavy-light meson chiral perturbation theory (HMxPT) to extrapolate Nf = 2 + 1 
lattice QCD results for 5-meson quantities to the physical quark masses. This follows the 
approach taken by the RBC and UKQCD Collaborations in the light pseudoscalar meson 



sector in Ref. 



34l |. and differs from the calculations of HPQCD and Fermilab/MILC, both 



of whom use SU{3) HMxPT for their chiral and continuum extrapolations ll5|. The 
use of SU{2) xPT is based on on the fact that the strange quark is much heavier than 
the up and down quarks, and can therefore be integrated out of the chiral effective theory. 
Because lattice QCD simulations at the physical strange-quark mass are possible via tuning, 
interpolation, or reweighting, SU (3) xPT is generally not needed to extrapolate the strange 
quark to its physical value. When the masses of the light valence and sea quarks are 
sufficiently small such that SU{2) chiral perturbation theory is applicable, SU{2) xPT for 
light pseudoscalar meson masses and decay constants converges more rapidly than SU{3) 



xPT |34l-l37j. Although we do not have enough data or sufficiently light quark masses to 
perform a thorough comparison of SU{2) and SU{3) HMxPT in this work, we believe that 
the use of SU (2) HMxPT provides a promising alternative to SU (3) and warrants further 
study when better data is available. 

In this work we compute both the ratio of decay constants, /s^/Zs^, and the ratio of 
AB = 2 matrix elements, ^. We focus on the S'f/(3)-breaking ratios because both the 
statistical and systematic errors are smaller and under better control than for the individual 
decay constants and mixing matrix elements. Our results have large total uncertainties 
compared to those of HPQCD and Fermilab/MILC. Within errors, however, our results for 
the SU (3)-breaking ratios are consistent with the values presented in the literature and we 
expect to improve upon them and present values for the individual decay constants and 
matrix elements in a future work. 

This paper is organized as follows. First, in Section Ull we present the actions and param- 
eters used in our lattice simulations. Next, in Sec. IIIIl we briefly discuss the perturbative 
matching of the heavy-light current and the four-fermion operators; the details of the lattice 



perturbation theory calculation will be presented in another publication 



We compute 
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2-point and 3-point lattice correlation functions and extract the decay constants and mixing 
matrix elements in Sec. llVt and extrapolate these results to the physical light-quark masses 
using a phenomenologically-motivated function based on next-to-leading order (NLO) SU (2) 
heavy-meson chiral perturbation theory in Sec. |Vl In Section I VI I we estimate the contribu- 
tions of the various systematic uncertainties to and ^, discussing each item in the 
error budget separately. We present our final results and conclude in Section IVIII 

This paper also contains four appendices. Appendix [X] specifies the SU (3) projection 
methods that are used in our APE and HYP-smeared gauge links. In App. [B] some details 
of the perturbative formulae used to match from continuum QCD to HQET as well as for 
the HQET running are presented. Appendix O discusses the large ground state degeneracy 
present in HQET and how it can be exploited to compute i?-meson mixing matrix elements 
using localized sources and sinks. The SU{3) and SU{2) NLO HM^PT expressions for the 
5-meson decay constants and mixing matrix elements relevant for Nf = 2 + 1 domain-wall 
lattice simulations are provided in App. [Dl some of these results have not been presented 
previously in the literature. 



II. LATTICE ACTIONS AND PARAMETERS 



In this section we briefly describe our numerical lattice simulations. We use the un- 
quenched lattices generated by the RBC and UKQCD Collaborations which include the ef- 
fects of 2+1 dynamical flavors of domain-wall quarks [l^. We calculate the decay constants 
and matrix elements on configurations with a lattice spacing of = 1.729(28) GeV 34| 
and an approximate spatial volume of ^ {1.8 fm)^. For each ensemble, the masses of the 
up and down sea quarks are degenerate and the mass of the strange sea quark is slightly 
larger than its physical value. In order to distinguish the dynamical quark masses used in 
our simulations from the physical u, d, and s-quark masses, we denote the lighter sea quark 
mass by mi and the heavier sea quark mass by rrih. Our lightest pion mass is approximately 
430 MeV. Table [T] summarizes the parameters of the dynamical domain- wall ensembles used 
in our analyses. 

In Section III Al we present the domain- wall fermion action used for both the valence and 
sea light u, d, and s quarks. Next, in Sec. Ill Bl we show the Iwasaki gauge action used for 
the gluon fields. Finally, we discuss the static action used for the heavy b quarks in Sec. Ill CI 
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TABLE I: Available 2+1 flavor domain-wall ensembles [l9|. The columns from left to right are the 
approximate lattice spacing in fm, the bare light and strange quark masses in the sea sector, the 
dimensions of the lattice in lattice units, the pion mass in the sea sector, the dimensionless factor 
m^L, and the residual quark mass in the chiral limit. 

a(fm) ami/anih Volume m7r(MeV) m^L amres 

« 0.11 0.01/0.04 16^x 32 430 3.9 0.00315 

« 0.11 0.02/0.04 16^ X 32 560 5.2 0.00315 

« 0.11 0.03/0.04 16^ X 32 670 6.3 0.00315 



A. Light-quark action 



We use the five-dimensional domain-wall fermion action 2l|, |22| for the light u, d, and s 
quarks in both the valence and sea sectors: 

'L,-l 



5l 



DW 



D\x,y) 



\s,s'=0 x,y X 

5^ ^ [(1 - l,)U,{x)6,+j:,y + (1 + 7^)f/;(2/)4-M,y] 



x] 



(5) 

(6) 
(7) 



{s = 0) 

(0 < s < - 1; 

(s = L, - 1) 



(8) 



where V's(^) is a 5-d Wilson-type fermion field. The fifth dimension extends from to 
Lg — 1 and is labeled by s. The domain- wall height is set to M5 = 1.8 in our simulations. 
The projectors Pl = (1 — 75)/2 and Pr = (1 + 75)/2 select left- and right-handed spinor 
components, respectively. The physical four- dimensional quark field q{x) is constructed from 
the five- dimensional field ips{.x) at s = and Lg — 1: 

qix) = PlM^) + Pr^l.-i{x), (9) 

qix) = %ix)PR + ^l._^{x)Pl. (10) 

In the limit — > 00, the left-handed and right-handed modes decouple and exact chiral 
symmetry is recovered. In practice, however, Ls is large but finite in numerical lattice simula- 



tions. This leads to a small amount of chiral symmetry breaking which can be parameterized 
in terms of an additive shift to the bare domain-wall quark mass. At the value Ls = 16 



used in our simulations we obtain a residual quark mass of arrij-cs = 0.00315(2) [3J]. Because 
mixing between heavy-light four-fermion operators of different chiralities is proportional to 
the value of arrives, this indicates that the size of errors from mixing with wrong-chirality 
operators is negligible. 

B. Gluon action 



We use the Iwasaki gauge action for the gluons |39|: 



'S'gaugc 



^(1 - 8ci) J2 I^eTr [Up] + ci J] ReTrff/^] j , (11) 



where (3 = G/g^ and go is the bare lattice coupling. Up is the path-ordered product of 
gauge links around the 1x1 plaquette P and Up is the path-ordered product of gauge 
links around the 1x2 rectangle R. The constant ci is set to —0.331 in the Iwasaki action 



41| for the quenched 



and we use (3 = 2.13 in our simulations. As was shown in Refs. |40L 
approximation, the use of the Iwasaki action in combination with domain- wall valence quarks 
leads to improved chiral symmetry and a smaller residual quark mass than for the Wilson 



gauge action 



42| . In the case of A'^^ = 2+1 dynamical domain- wall simulations, the use of 



the Iwasaki action also allows frequent tunneling between topological sectors 



C. Heavy-quark action 

We use an improved static action for the 6-quarks in the 2-point and 3-point correlation 
functions needed to compute the decay constants and matrix elements. We build upon the 



origma. 
Hill 



4|: 



lattice formulation of the static effective action that was constructed by Eichten and 



Astatic = a'Y, {H^)[^^,y - Ul,{y)S,.ojP+h{y) - h{xM,,y - [/o(x)<5,^o,J^-^(?/)) , (12) 

x,y 

where h{x) is the static quark field at site x, Uo{x) is the gauge link in the temporal direction 
and denotes the unit vector along the temporal direction. The projectors P± = |(1 ± 70) 
select the parity even and odd components of h{x), which we denote by h^^\x) and h^~\x). 



respectively. The components of the static-quark field h^"^^ in the Eichten-Hill action satisfy 
the relation 7o/i''^^ = i/i'-^'*. Furthermore, the static quark propagator H{x, y) = {h{x)h{y)) 
can be expressed as the product of gauge links: 



H{x,y) = H'^{x,y) + H {x,y), 

y) = ^9{t, - ty)5,,g [Ul{x - 0) . . . f/J(y + Q)Ul{y) 
H^~\x, y) = -\e{ty - Q6s,^ \Uoix)Uoix + 0) . . . U^iy - 0) 



(13) 
(14) 

(15) 



These properties make the Eichten-Hill formulation computationally simple. In practice, 
however, numerical simulations with this action are quite noisy. Therefore we use instead 
smeared (or "fat") link actions: 

5fat = a'J2 {hi^)[Ky - vl{y)5,_o,y]P+hiy) - h{x)[5.,y - Vo{x)5,^^,jP^h{y)) , (16) 
x,y 



27 



28 



where the new gauge link V is obtained from the thin link U by either APE blocking 
or hypercubic (HYP) blocking j29|. The heavy quark propagator for the improved action is 
given by 



T7t, 



T7t, 



H\:h^^ y) = --,d{ty - U)6s,y- Vo{x)Vo{x + 0) . . . Vo{y- 0) 



(17) 
(18) 



The replacement of the simple gauge link Uq by a smeared link significantly improves the 



signal-to-noise ratio 



26|. 



We construct the fattened APE link by adding a weighted sum of the staples to the 
original thin link, and restrict the smearing to links in the temporal direction along which 
the heavy quark propagates. We use the APE parameter a = 1, for which the smeared link 



is given by 



Vo{x) = Projgu(3)[K)(x)], 



(19) 



1 

Voix) = -Y,{M^W^ + ^)Ulix+d) + Ulix-u)Uo{x-u)U,{x+d-u)y (20) 

We build the HYP link from three iterative steps of APE smearing which are restricted to 
the hypercube around the original link. For links in the temporal direction the construction 
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is as follows: 



Proj 



SU(3) 



1 - «l)f/o(x) + ^ ^ K;0(x)K);.(x + i>)V;yX + 0) 



±3 



Proj 



SU(3) 



(1 - a2)U,{x) + ^Y. 



v=±l 
±3 



p=±0 



V;.;i.p(a;) = Proj 



SU(3) 



[l-a,)U,{x) + '^ J2 Ur,{x)U,{x + ri)Ul{x + fl) 



±3 



r;=±0 



(21) 



(22) 



(23) 



We use the HYP smearing parameters (ai, 0^2, 0^3) = (1.0, 1.0, 0.5), sometimes referred to as 
HYP2. These were shown to approximately minimize the noise-to-signal ratio in Ref. 25 1. 

In Eqs. (fT9i)-(l23l). Projsu(3)(V^) indicates the projection of the link V onto an SU{3) 
matrix. This projection reduces the statistical noise and thus enhances the smearing effect 

without increasing the level of smearing. In some cases the S'?7(3)j3rojection also suppresses 

I, 

quantum corrections to lattice operators in perturbation theory [45|. The projection of the 
smeared link onto SU{3) is not unique, and we use two different schemes: for the case of 
APE smearing, Eq. ( fT9l) . we project by the unit circle method based on polar decomposition 
461], while for the case of HYP smearing, Eqs. (l2T]) -(!23l). we obtain the projected matrix 
by an iterative procedure seeking the SU(3) matrix f/max that maximizes Re Tr (f/max^^), 



where V is the HYP smeared link matrix 



47| . We describe the details of the two schemes 



and show their equivalence in the weak coupling limit in Appendix |Al 



III. PERTURBATIVE MATCHING OF HEAVY-LIGHT CURRENT AND FOUR- 
FERMION OPERATORS 

In order to renormalize the heavy-light axial current and AB = 2 four-fermion operator, 
we adopt a two-step matching procedure. In the first step, we match the QCD operators 
renormalized in the MS scheme using naive dimensional regularization (NDR) at a scale 
fib to continuum static effective theory operators renormalized at a scale /i. This step is 
described in Sec. IIII A[ In this paper, we choose fib to be the 6-quark mass nib and fi to be 
the inverse lattice spacing a~^. In the second step, we match the continuum static effective 
theory operators to the lattice ones. This step is described in Sec. IIII B[ We combine the 
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results of the two steps and present the results for the complete matching coefficients in 
Table [ml 



A. Continuum matching 

The QCD operators considered in this paper are the axial vector current 

A?"""" = 670759, (24) 
and the AB = 2 four-quark operator 

0?^°(/i,) = [h,{l - l,)q][h,{l - 75)g]. (25) 
These are related to the continuum HQET operators by 

= C^(/i)<Q'^^(/i) + 0(AQCD/m,), (26) 

0?^°(/i,) = /i)0^«^^(/i) + Z2(/i,, /i)0?Q''^(/i) + 0(AQCD/mb), (27) 

where 

A^Q^T = 7^70759, (28) 

O^Q^^ = - 75)g]fc(l - 75)g], (29) 

o^Q'^^ = Ri-75)g]Ri-75)g]. (so) 

Note that, because A^*"^ is a conserved current, it does not depend upon the renormalization 
scale fib- The coefficients Ca and Z = {Zi, Z2) in Eqs. (126|) and (127]) are products of three 
factors: 

= CA{mk) ■ U^^\mb,mc) ■ uf^{m^,^), (31) 
Z(/ib,/i) = Z{iJ,b,mb) ■ U]^\mi„mc) ■ U^^\mc,^), (32) 

where Ca and Z = (Zi, Z2) are the matching coefficients from HQET to QCD. The factors 
U^J^^\n', /i) and the 2x2 matrix U^^\n', fx) account for the renormalization group running 
between scales: We first match the continuum QCD operators onto continuum HQET oper- 
ators at one-loop; this occurs clt cL SCclIg fib = TUb in the MS(NDR) scheme. We then run the 
matching coefficients in four-flavor continuum HQET from rrib to rric at two-loops. Finally, 
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we run the coefficients at two-loops in three-flavor continuum HQET from rric to the scale 
/i = 1/a, where the matching to three-flavor lattice HQET is done. 

The matching factors and anomalous dimensions needed to compute these coefficients are 



given in Refs. 



30 



48 



-|53| for the four-quark 



49| for the heavy-light current and in Refs. (50- 
operator. For completeness, we present them in Appendix |Bl Here we simply quote the 
results for the matching coefficients, which already contain some terms of 0{a^): 



CAifi) 



asirrib) 


6 

25 











1 - - 



as{mb) 



3 Air 



Zi(mb,/i) 



i + 




(33) 



1 + j(3)^("^c) - 



4:71 



Z2(m6,/x) = -8 



as[mb) 
47r 



as{mb) 


25 











(34) 
(35) 



with the parameters Ja and Jn given by 



J 



(3) 



-0.7545, 



jJJ) = -1.6980, 



t(4) 

-'a 

•^11 



-0.9098, 
-1.8637. 



(36) 
(37) 



Note that the leading-order mixing between 0^^^ and Og'^^^ is of 0{as). To determine 
the coupling constant as at different scales, which is required to obtain these results, we flx 
the value of ag at the Z-boson rnass to the PDG value aginiz = 91.jL876 GeV) = 0.1176 j^. 



54 



Using four-loop running 
1.27 GeV jfil) = 0.3819, and a 



55| . we obtain a.^mh = 4.20 GeV 



(3)/ -1 

'a 



0: 



0.2228, ajrur. 



1.729 GeV [3^) = 0.3141. Hence we flnd 



CA{a^^) = 1.0459, Zi{mb,a^^) = 0.9100, Z2{mb, 



-0.1502 



(38) 
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for the matching coefficients in Eqs. (126|) and (127|) that relate the operators in continuum 
QCD to those in the continuum static effective theory. 



B. Static effective theory matching 

We next discuss the lattice-to-continuum operator matching in the static effective theory. 
Although both the domain- wall and static quark actions are on-shell 0{a) improved, the 



domain- wall-static vertices still receive 0{a) corrections. This was shown in Ref. [561] for the 
case of clover light quarks and static heavy quarks even when the parameter in the clover 
action r — )■ and chiral symmetry is restored in the light-quark sector. In this work, we 
remove 0{a) lattice discretization errors from the axial vector current in the light-quark 
chiral limit at one-loop in lattice perturbation theory by adding dimension-four operators 
containing derivatives; we refer to these errors as 0{pa) to distinguish them from 0{mga) 
errors that vanish in the light-quark chiral limit. The improved operators are determined 
by requiring that the on-shell quark scattering amplitudes in the 5-meson rest frame agree 
in the lattice and continuum theories through 0{aspa), where p is the momenta of the 
light quarks in the 5-meson and is typically of (9(Aqcd)- The details of this calculation 



are presented in Refs. 38|]. We neglect 0{mga) errors because these are estimated to be 



small 381], but account for them when estimating the systematic errors in Sec. IVI[ In this 
work we also neglect 0{aspa) corrections to the four-quark operator since we estimate the 
size of the resulting contribution to the SU{3) breaking ratio C, would be much smaller 
than our current statistical precision, but we again account for them in the systematic 
error budget. The one-loop 0{aspa) corrections for the four-quark operator are available in 



Ref. 



38 



however, for use in future simulations. 



For domain-wal 



presented in Refs. [32 



it quarks, the one-loop lattice perturbation theory calculations are 



38 



57(1 . The lattice-to-continuum matching has the form: 
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with the constant 

<^ = l-M5 + 4(l-no), (42) 

where the mean-field hnk uq = P^^^ is obtained from the expectation value of the plaquette 
P. In Eqs. dSHD-dH]) we match A^^^^ and O"^^^ at one-loop, but match O^^^'^ at tree-level. 
This is sufficient because the leading-order mixing between the continuum QCD operator 
Cf^^ and the continuum HQET operator Og'^^^ is already of 0{as), and has no tree-level 
component. The lattice operators have the same form as in the continuum static effective 
theory 

K^' = hlol5q, (43) 
Of = [%(1 - 75)g][/^7/.(l - 75)g], (44) 
Of = [h{l-^,)q][h{l-^,)q]. (45) 

The 0{pa) derivative operator in the equation for the axial current is given by 

A'Ho = do (hl.q) , (46) 

where we have simplified the expression using the equations-of-motion. The domain-wall 



specific renorma 



ization factor that enters the above equations was calculated pertur- 



batively in Ref. 58 1. 

The superscript "MF" denotes mean- field improvement [59|] , in which we modify the bare 
lattice coupling using the value of the mean-field link. Use of this "boosted" coupling as 
the new expansion parameter improves the convergence of lattice perturbation theory. The 
mean-field improved coupling has several definitions which differ only at higher-order in 
perturbation theory than we consider. These differences enter our estimate of the systematic 
error in Section I VII Our choice for obtaining the mean-field improved coupling from the 
bare lattice coupling is 

1 P 

= ^ + + + ^/^/' ('^^) 

where Nj = 3 is the number of dynamical fiavors. We use the plaquette value averaged 
over different light sea quark mass, P = 0.5881, because the difference in the value of P 
for different light sea quark masses is less than 0.05% for the ensembles used in this work. 
The constants dg = 0.1053 and Cp = 0.1401 were calculated for the Iwasaki gauge action in 
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TABLE II: Lattice-to-continuum operator matching factors in the static effective theory for the 
choice of coupling . 



smearing Zf^ c^^ Zf^ Zf^ 
APE 0.9090 0.0653 0.8225 0.9642 
HYP 0.9382 0.1204 0.8909 0.9642 



Ref. 58[ and df = —0.001465 was obtained for = 1.303 from a linear interpolation of 



60|. On our ensembles, the mean- field improved coupling is 



the results in Table II of Ref. 
^MF ^ 0.1769. 

Given the value of the plaquette in our simulations, uq = 0.8757 and = —0.3029. 
After setting yU = a~\ the coefficients appearing in Eqs. fl39l) . fj40|) and fHTl) are: 

= 1 + ^3 X 5.250, (48) 

c^MF. I -1.584 APE 
Zr = 1 + X I (49) 



47r 3 



0.077 HYP 



a^F4 3.480 APE 
I 6.412 HYP 

MF -4.462 APE 
I 1.076 HYP 

where the values are given for both the APE and HYP link-smearings used in this work. 
These results can be combined to determine the overall multiplicative renormalization factors 
denoted by Z in Eqs. ( I39l) - (HT1) . which we present for completeness in Table HTl 



C. Complete matching coefficients 

In order to match the lattice HQET operators at scale /i = directly onto the desired 
continuum QCD operators at fii, = rrih, we must combine the coefficients obtained in the 
two steps. We define the complete lattice HQET-to-continuum QCD matching coefficients 



15 



TABLE III: Perturbative matching factors for the decay constants and mixing matrix elements 
evaluated for APE smeared and HYP smeared static-quark gauge links for the choices of the 
strong coupling constant . 



smearing 


; Z$ ZvA 


Zsp 


APE 


0.9507 0.7485 


-0.1448 


HYP 


0.9813 0.8108 


-0.1448 



as 

Z^{a-') = CA{a-') ■ ZTi<^~\ a-^), (52) 

a"') = Ziifi,, a~') ■ ^(a"\ a^^), (53) 

Zspifib, a-') = Z^ifib, a"') ■ «^'), (54) 

and present their values for our choice of simulation parameters in Table lllli These will 
be used in the following section to extract the physical decay constants and AB = 2 four- 
fermion matrix elements via the relations 

= Z^{a-') «*(a-i) + cTaA'g%a-')) , (55) 
0?^°(/i,) = ZvAifib, a-') Of (a-i) + Zsp{f^b, a"') Of (a"^). (56) 

In practice, the renormalization factor Z^ cancels in the ratio fss/ fB^i and only the 
quantity Zsp/Zva enters the ratio ^. Therefore we do not need Zq, (or Zva and Zsp by 
themselves) for our current analysis of the S'f/(3)-breaking ratios. We present all three 
matching coefficients for completeness, however, because they will be necessary for calcu- 
lating the individual decay constants and four-fermion operator-mixing matrix elements in 
future work. 

IV. LATTICE CALCULATION OF 5C/(3) BREAKING RATIOS 

In this section we calculate the ratios of the 5-meson decay constants and mixing matrix 
elements at unphysical values of the light and strange quark masses. On each sea quark 
ensemble, we compute the necessary 2-point and 3-point correlation functions at two values 
of the valence quark mass: the unitary point rrix = mi and a point tuned to the physical 
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TABLE IV: Parameters used in our simulations. The columns from left to right are the light and 
(approximately) strange sea quark masses, the light and strange valence quark masses, and the 
number of configurations analyzed for both of our setups using APE or HYP link smearing in the 
static quark gauge links. 









# configs. 


ami/anih 


amx 


APE HYP 


0.01/0.04 


0.01, 


0.0359 


298 300 


0.02/0.04 


0.02, 


0.0359 


298 300 


0.03/0.04 


0.03, 


0.0359 


298 300 



strange quark mass anis = 0.0359 [19|. We also use two different link smearings (APE 
and HYP) to improve the static heavy quark action in order to help estimate discretization 
effects. Table IIVI presents the parameters chosen for our matrix element computations. In 
the first subsection we calculate the ratio of 5-meson decay constants and in the second we 
calculate the ratio of AB = 2 mixing matrix elements. 



A. Calculation of the ratio of i?-meson decay constants 

In QCD the decay constant /b^ for the iJ^-meson is defined by the vacuum-to-meson 
matrix element 

{0\h,l5q\Bgip))=tfB,Pt.. (57) 

Because the decay constant /b, behaves as 1/ y/ruBq in the Umit of large S^-meson mass, we 
calculate the combined decay amplitude 

= fB,,/m^, (58) 

where is the physical mass of the S^-meson. We determine the quantity ^Bg by com- 
puting two-point correlation functions of the static-light axial current j\'-^'>^^^^ = h^'^^^'j^q} 

^ For the case of static 6-quarks, we can relate the i?-nieson interpolating operator ''75'? to the axial 
current operator in the temporal direction using the relation /i(^-'(a;)7o = ±/i(^-'(a;), and thereby express 
all correlation functions entirely in terms of the axial current. 
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In practice, we use Coulomb gauge-fixed wall sources for the 6-quark to calculate the 
local-wall (LW) and wall-wall (W^iy) correlation functions 



C^^(t,to) = a'Y.^O\A^(x,t)A^ito)^\0), (59) 

x&V 

C^^(t,to) = (OK(tX(to)^lO), (60) 

with the local (L) and wall-source (W) axial currents given by 

Ao(f,t) = 7i^"^\f,t)7o75g(x,t) +7i^ ^(f,t)7o75g(x,t), (61) 
A^it) = a^J^J^ (7^^"'\y,t)7o75g(^"',t) + /^^"\y,t)7o75g(^"',t)) • (62) 

From the ratio of to C^'^ we obtain the combined decay amplitude 

Uat /2 \C^^{t,to)\ 




= lim W — ' ^ ' (63) 
where we determine the unphysical S-meson rest mass m^^ via 

A derivation of Eq. (!63|) is presented in Appendix Finally we compute the renormalized 



decay amplitude 



= [l + c7 sinh (am^ J J $1^;, (65) 

using the perturbative matching factors given in Tables HTl and lllll The contribution propor- 
tional to c^^ improves the heavy-light axial current operator through 0{aspa), where the 
sinh arises from the symmetric derivative in the 0{pa) operator. The overall multiplicative 
factor Z$ is needed to obtain the combined decay amplitude in the continuum. 

The computation of the statistical errors throughout this paper follows the prescrip- 



tion for numerically computing the autocorrelation function as proposed in reference |61 |. 
The autocorrelation function quantifies the degree of correlation between two measurements 
made at different trajectories, and depends upon the observable of interest. By summing 
the autocorrelation function over the separation between measurements, one obtains the 
integrated autocorrelation time. We obtain a better estimate of the true statistical error by 
inflating the variance of the measured Monte Carlo data using the integrated autocorrela- 
tion time. In many instances we must compute the errors in a quantity which itself depends 
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on several lattice correlators; we refer to this as a derived quantity and refer to the lattice 
correlators as primary observables. When calculating the errors in a derived quantity, we 
account for the correlations between primary observables by using the functional dependence 
of the derived quantity on the primary observables. As a cross-check of the statistical error 
estimate, we compared the results obtained with this approach with those obtained using a 
single-elimination jackknife procedure; we find that both the central values and statistical 
errors are consistent between the two error estimation methods. 

For example, Fig.[T]shows the determinations of m^^ and on the ami = 0.02 ensemble 
for the APE (upper plots) and HYP data sets (lower plots). The central value and statistical 
error of each data point in Fig. [T] are computed as functions of the primary observables 
and . Then the value of the plateau and its error are computed using a function which 
averages the values of fn*^^ (or $^") on time-slices 12, 13, 14 and 15 because we do not observe 
excited-state contamination in this region. In order to reduce the size of the statistical errors, 
we average the correlators beginning at two time sources. We achieve the averaging of our 
two sources by replacing e.g. C^^{t, 0) with [C^^{t, 0) + C^^(20a - t, 20a)] /2 in Eqs. ([63])- 
(165|) . For the case of the HYP-smeared data, these are located at t/a = and t/a = 20, 
whereas for the APE-smeared data the second source is located either at t/a = 20, 21, or 
24. 

Finally, we compute the ratio $^^/$^° on each ensemble as a function of and 
which themselves depend on the corresponding primary observables and C^^-^ this 
is shown in Fig. O We obtain the plateau for the ratio from time-slices 12-15, where we 
do not observe excited-state contamination in the numerator or the denominator 
Table IVl presents the values of ^^"/'l''^" on the three sea-quark ensembles. Despite the use 
of two time sources, the statistical errors in the ratio are as large as 7.5% in case of the data 
using APE smearing. 

B. Calculation of the ratio of i?-meson mixing matrix elements 

The Bg — Bq mixing parameter in continuum QCD is defined in terms of the matrix 
element of the AB = 2 four-fermion operator via Eqs. ([2]) and ([3]). Because the matrix 
element Aig behaves as in the limit of large i?g-meson mass, we calculate the desired 
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FIG. 1: Determination of (left) and <I>^° (right) using the average of two time sources on 
the ami = 0.02 ensemble. The upper plots show the APE data, while the lower plots show the 
HYP-smeared data. For each panel, the shaded band corresponds to the plateau extracted from 
averaging the data over four consecutive time slices. Errors shown are statistical only. 



matrix element divided by the i?g-meson mass: 



(66) 



For the determination of Mb^ , we use different spatial sources for the two choices of link 
smearing: in the case of APE-smeared links we use box sources of size 8'^, while for HYP 
smearing we use wall sources. The local-wall two-point functions are already defined in 
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FIG. 2: Determination of $^°/$^° using the average of two time sources on the three sea quark 
ensembles. The blue (triangle) points denote the APE data, while the red (square) points denote 
the HYP-smeared data. The shaded (hatched) band corresponds to the plateau extracted from 
averaging the APE (HYP) data over four consecutive time slices. Errors shown are statistical only. 



Eq. (!59|) . while the box-box two-point functions are given by 



C^^(t,to) = (0|A-(t)A-(to)^|0) 



(67) 



with 



(68) 



x,yeAV 

where the superscript B denotes a box source in the region A^. We also compute the 
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TABLE V: The renormalized decay amplitude ratio <^^^/<l>^° for both the APE- and HYP-smeared 
data sets. Errors shown are statistical only. 







/$-n 


ami 1 CL^'t-h 


APE 


HYP 


0.01/0.04 


1.07(8) 


1.10(4) 


0.02/0.04 


1.06(5) 


1.02(3) 


0.03/0.04 


0.97(4) 


0.99(2) 



three-point correlation functions 

CS.(t;,t,to) = a3 5^(0|4^)'(t;)tor(f,t)4"^'(to)^|0), (69) 

where A^^''^{t) is either the box-source axial current A^f^''^{t) or the wall-source axial current 
Al^^^{t). The operator Oi can be either the sum of the squared vector plus squared axial 
vector current or the squared scalar plus squared pseudoscalar current, respectively,^ 

O'vWaa = 2 {h^^\'q) (h^~\,q) + 2 (h^^\^l,q) (h^'^^l^q) , (70) 
Ofs+PP = 2 (h^^\) (h^-\) + 2 (h^^\,q) (h^~\,q) • (71) 

Although the SS + PP operator does not contribute to Bg — Bq mixing in continuum QCD, 
its counterpart in HQET is introduced through the QCD — t- HQET matching as shown in 
Eq. ([27D. 

Because of the different spatial wavefunctions used for the APE- and HYP-smeared data, 
we extract the Bg — Bq matrix element in different ways for the two data sets. With the 
APE-smeared box source data, we can compute the matrix element directly from the ratio 



of correlators 



57| 



Mo- = lim 2 ^ ' 72 



We neglect the odd-parity parts of O^* and Og*, Eqs. ([^^ and (P5|). when computing the lattice three- 
point correlation functions because only the parity-conserving components contribute to the desired matrix 
element. 
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For the HYP-smeared wall-source data, however, we must first extract the bag parameter 
Bb 

3 3 ^^^fl^l^^ 

= t,l%u, 8^ C^^(t,t;)C^^(t,to)' ^^^^ 
then we obtain the matrix element via 

Mo, = ^-Bo,{n:f- (74) 

The derivation Eq. ( 1721) is more complex than that of either Eq. (l63l) or Eq. ( 1731) because 
the box sources in the amplitudes appearing in both the numerator and denominator are 
not translationally invariant and will create 5-meson states carrying a variety of spatial 
momenta. As explained in greater detail in App. [Cl these states are degenerate with the 
lowest energy 5-meson state and so cannot be suppressed by simply going to large Euclidean 
time separations. Instead, Eq. ( !72|) is derived in App. Ousing the local conservation of heavy- 
quark number, a property special to the static approximation. Finally, for both the APE 
and HYP-smeared data, we compute the renormalized matrix element through 0{as)i 

M-'^ = Zyj^M'^'y^^^ + ZspMl%pp, (75) 

using the perturbative matching factors given in Table IIIII 

We compute the ratio of mixing matrix elements in the same manner as we compute 
the ratio of decay constants in the previous subsection. For example. Fig. [3] shows the 
determination of the numerator M^°J^ and the denominator M'^^^ on the ami = 0.02 ensemble 
for the APE (left-hand plot) and HYP data sets (right-hand plot). We compute the values of 
the plateaux using time-slices 8-12 because we do not observe excited-state contamination 



in this region. Table |VT] shows the results for the S'f/(3)-breaking ratios ■^M^^/M'^^ on the 
three sea quark ensembles; the corresponding plateau plots are shown in Fig. HI 



V. CHIRAL EXTRAPOLATION 



We extrapolate our results for the ratio of the decay constants and for the ratio of mixing 
matrix elements to the physical point using a phenomenologically-motivated function based 
on next-to-leading order partially quenched SU(2) heavy-light meson chiral perturbation 
theory. In SU{2) xPT, the pesudoscalar mesons containing strange quarks (i.e. kaons and 
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FIG. 3: Determination of M^'^ and M^^^ on the ami = 0.02 ensemble. The left-hand plot shows 
the APE data, while the right-hand plot shows the HYP-smeared data. For each panel, the shaded 
band corresponds to the plateau extracted from averaging the data over four consecutive time 
slices. Errors shown are statistical only. 

TABLE VI: The renormahzed 5C/(3)-breaking ratio ^M'gf/M'g^ for both the APE and HYP- 
smeared data sets. Errors shown are statistical only. 





jMf;/M- 


ami/amh 


APE 


HYP 


0.01/0.04 


1.050(78) 


1.110(49) 


0.02/0.04 


1.038(40) 


1.006(38) 


0.03/0.04 


0.992(27) 


0.987(17) 
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r/'s) are integrated out of the theory. Thus SU{2) xPT does not require an expansion in 
the strange quark mass about the chiral Umit, and the SU{2) xPT expansion parameter 
in isospin-symmetric simulations is rrii/A^, where m; is the light up-down sea quark mass 
and is a typical hadronic scale. This improves the convergence of the chiral expansion 
relative to SU{3) xPT, as long as mi is sufficiently light that corrections of 0{mi/ms) 
are small. Studies by the RBC and UKQCD Collaborations, the PACS-CS Collaboration, 
and the MILC Collaboration confirm this picture and show that, for light pseudoscalar 
meson masses and decay constants, SU (2) xPT within its applicable region converges more 
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FIG. 4: Determination of ^ JVP^"^ / K'P^^ using the average of two time sources on the three sea 
quark ensembles. The blue (triangle) points denote the APE data, while the red (square) points 
denote the HYP-smeared data. The shaded (hatched) band corresponds to the plateau extracted 
from averaging the APE (HYP) data over four consecutive time slices. Errors shown are statistical 
only. 



quickly than 5'f/(3) xPT |34l-l37l|. For the case of S'[/(3)-breaking ratios such as fsJlBa 
and ^, however, SU{3) HM^FT has the advantage that the chiral extrapolation formulae 
manifestly preserve the fact that the ratios must be equal to one in the limit mi — > m^. 
Within the framework of SU(2) HM^FT, this fact must be introduced in a more ad hoc 
manner such as by matching the SU (2) HMxFT expression at small quark masses onto an 
analytic form at large quark masses that becomes one when m/ — )■ m^. We therefore plan 
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to perform a study comparing the use of SU{2) versus SU{3) HM^PT for the extrapolation 
of i?-meson decay constants and mixing matrix elements to the physical quark masses in 
a future analysis when we have lighter data with smaller statistical errors in order to see 
which procedure leads to a more accurate determination of these quantities. 

Although we know that HMxPT is the correct low-energy effective description of the 
lattice theory when the simulated quark masses are sufficiently light, we do not know a priori 
at what mass the range of validity of HMxPT ends. Studies of the light pseudoscalar meson 
sector on the RBC/UKQCD domain- wall ensembles show that NLO xPT does not describe 
the numerical data for the masses and decay constants when the pion masses are above 



about 420 MeV [3J] . Since the lightest pion mass in our analysis is approximately 430 MeV, 
this suggests that most of our data may be too heavy for NLO HM^PT to apply and that 
the inclusion of NNLO terms may be necessary. Unfortunately, because we only have three 
data points for each of the S't/(3)-breaking ratios, we do not have enough data points to 
reliably constrain the values of the higher-order terms (there are two free parameters at 
NNLO in the SU{3) HM^PT expressions, and even more at NNLO in SU(2) HM^PT). 
When the masses of the light pseudoscalar mesons in the chiral logarithms are sufficiently 
heavy, however, the logarithms can be well-approximated by polynomials. We therefore 
choose to fit the data to a linear fit function that is constrained to be equal to one in the 
SU (3) limit (m; — )■ m^) and then match this result onto the NLO SU (2) HM^PT expression 
at the location of our lightest data point {ami = 0.01), which we believe is sufficiently light 
that NLO SU{2) HMxPT should apply. Although the behavior of the data is unlikely 
to be strictly linear in the heavy-mass region, the statistical errors in our data points are 
sufficiently large (as great as ~ 8%) that we can successfully perform a linear fit and obtain 
a good x^/dof without the addition of higher-order polynomial terms. ^ Furthermore, we 
cannot reliably determine the size of a quadratic term if we include one in the fit. Given 
our poor statistical errors, however, we cannot exclude the possibility of other fit functions, 
and we use alternate fit forms as one way to estimate the chiral extrapolation error. We 
also vary the location of the SU (2) HMxPT matching point and the parameters that enter 
the SU (2) HM^PT expressions, and consider matching onto SU (3) HMxPT. All of these 



^ Because the data points for the SU (3)-breaking ratios in Tables IVl and I VII were generated on three different 
sea-quark ensembles, they are statistically independent; thus the uncorrelated x^/dof correctly reflects 
the goodness-of-fit. 
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variations are discussed in greater detail when we estimate the systematic uncertainty due 
to the chiral extrapolation in Sec. IVI A[ 

In the case of the decay constants, we extrapolate the ratio 



while for the mixing matrix elements we consider 



Ml 



(77) 



From here on, we drop the superscript "ren" for simplicity because we only refer to the renor- 
malized quantities. The expressions for the heavy-light meson decay constant and mixing 
matrix element to NLO in the light-quark expansion, but zeroth order in l/rrib, are given in 
Appendix For completeness, we present formulae for both SU (3) and SU (2) HMxPT. 
The SU{3) fit functions depend on the valence and sea light-quark masses {ms,mi,mh) and 
the lattice spacing a. The SU{2) functions are obtained from the SU{3) expressions by 
taking the limit ( ^^^] <^ 1. Thus they only apply in the region in which the value of 
the average up-down quark mass is much smaller than the valence and sea strange quark 
masses. Furthermore, because the strange quark has been integrated out of the SU (2) the- 
ory, the SU (2) fit functions only depend upon the light-quark mass mi and the lattice spacing 
a. Therefore the expressions for the S'?7(3)-breaking ratios at NLO in SU(2) HMxPT are 
particularly simple: 



■SfB, n (lirff \1J \AIJ '(to/) 



where m| = 2B{mi + m.^cs) is the tree- level mass-squared of a pseudoscalar meson composed 
of two quarks with mass m/, and / are the leading-order low-energy constants of xPT, 
and the quark masses in the analytic terms are expressed in terms of dimensionless ratios to 
make the coefficients Ci and Di of 0{1). These functional forms are derived from Eqs. ( 1016^ - 
(ID20p by taking the ratio of the expressions for the valence quark y = s over the expressions 
for X = I. Because we are working at a single lattice spacing, the analytic terms proportional 
to are absorbed into the values of the leading-order coefficients and Rm- Note that 
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in the limit mi —¥ rUs the 5'[/(3)-breaking ratios are not constrained to unity, as would be 
the case in SU{3) HMxPT. (In fact, the point mi = does not even lie within the range 
of validity of SU(2) HMxPT and hence of Eqs. (1781) and (17^ .) This is because, once the 
strange quark has been integrated out of the SU (2) theory, the expressions no longer contain 
explicit strange-quark mass dependence. All of the effects of the strange quark are encoded 
in the values of the low-energy constants, which differ in the SU (2) and SU{3) theories. 

Although the coefficients of the chiral logarithms depend on the low-energy constants 
gB*B-Ki f 1 and S, once these are fixed as we now describe, there are only two free parameters 
each in Eqs. ( ITHj) and ( 1791) : the overall normalization and the coefficient of the analytic term 
proportional to m^. This allows us to smoothly match the SU(2) expressions onto the linear 
fit of the heavy data without ambiguity. In the chiral extrapolation we obtain our central 
value using gB*BTT = 0.516 for the B*-B-7r coupling, which comes from a two-flavor lattice 
determination in the static heavy quark limit by Ohki, Matsufuru, and Onogi 62] • We then 
vary the value of gB*Bn over a reasonable spread of values based on both lattice calculations 
and phenomenological fits to experimental data in order to estimate the systematic uncer- 
tainty, as described in further detail in Sec. IVI B[ Moreover we set the leading-order pseu- 
doscalar meson decay constant / to the experimental value of f.^ = 130.4±0.04±0.2 MeV j^. 
This is consistent to the order in xPT at which we are working since it only modifies higher- 
order NNLO terms. Studies by both the MILC and JLQCD Collaborations suggest that the 
use of a phy sical parameter in the chiral coupling (/ — t- /^) leads to improved convergence 



of xPT 63|, |64] . The scale in the chiral logarithms is fixed by setting A^^ = 1 GeV. For the 



low-energy constant B we use the value aB = 2.414(61) obtained from a NLO fit of the 
pseudoscalar meson masses 3J|. Finally, whenever the residual quark mass appears, we use 
its value in the chiral limit amres = 0.00315. 

The results of the chiral extrapolation are shown in Figs. |5] and |6l The blue triangles (red 
squares) show the data obtained using APE (HYP) link smearing and are plotted versus the 
light sea quark mass. We indicate the location of the physical strange quark mass by 
the black dot. The dashed vertical line marks the physical average u-d quark mass, which 
is the point at which we extract the physical values for ^Ss/^Ba ^^'^ a/ ^Bs/^b^ ^- The 
agreement between the two smearings is good. For the case of the APE data, the x^/dof 
for the fit of both S'f/(3)-breaking ratios is below one, indicating that the data are well- 
described by the linear fit function. For the HYP data, the x^/dof's are 1.8 for ^bJ^Bj, 
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and 2.0 for \JmBj'mB^ ^, respectively. These still correspond to confidence levels of greater 
than 10%, however, and are therefore consistent with the data. The error bands in Figs. [5] 
and [6] are fairly broad and hence the statistical uncertainty at the physical point is large, ~ 
4.1-7.6% for (^bJ^b^ and ~ 4.7-6.3% for ^JmBjraB^ i- This is due to the large statistical 
errors in our data points, as well as the fact that our lightest mass is still quite heavy, which 
forces us to extrapolate over a large mass range. We summarize the results of our preferred 
chiral extrapolation in Tab. IVIII and discuss the estimation of our systematic errors in the 
following section. 



APE: <i>Bj'iBi = 1.165(88), xVdof = 0-8 
HYP: '^bJ'^b, = 1.153(47), xVdof = 1.8 
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FIG. 5: Chiral extrapolation oi ^Bs/^Bi = yrnjiJ^nBi • fBs/ Ibi- The blue (triangle) points 
denote the APE data, while the red (square) points denote the HYP-smeared data. The color 
of the shaded (hatched) error band corresponds to those of the APE (HYP) data points. The 
dashed vertical line denotes the physical average u-d quark mass and the black dot denotes the 
physical strange quark mass, at which the S'C/(3)-breaking ratio must be one. The SU (2) HM^PT 
coefficients obtained from the fit are = 1.21(9) and C/ = 1-1(6) in the case of APE smearing 
and = 1.19(5) and Ci = 1.2(3) in the case of HYP smearing. Errors shown are statistical only. 
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APE: fnBjniB, ■ i = 1.142(72), x'/dof = 0.3 
HYP: ^jmBjniB^ ■ ( = 1.144(54), xVdof = 2.0 
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FIG. 6: Chiral extrapolation of ^^MbJMbi = \fm^Jra^^ (/Bs \J ^Bs I fsi a/ Bbi)- The blue (tri- 
angle) points denote the APE data, while the red (square) points denote the HYP-smeared data. 
The color of the shaded (hatched) error band corresponds to those of the APE (HYP) data points. 
The dashed vertical line denotes the physical average u-d quark mass and the black dot denotes the 
physical strange quark mass, at which the S'[/(3)-breaking ratio must be one. The SU (2) HM^PT 
coefficients obtained from the fit are Rm = 1-18(8) and Di = 2.6(1.1) in the case of APE smearing 
and Rm = 1-18(6) and Di = 2.6(8) in the case of HYP smearing- Errors shown are statistical only. 



TABLE VH: Results for ^Bs/^b^ and yrnsJrnB^ - at the physical point. Only statistical errors 
are shown. 





link smearing 




APE HYP 


^bJ^b, 


1.165(88) 1.153(47) 


^JniBjniB^ 


^ 1-142(72) 1-144(54) 



VI. ESTIMATION OF SYSTEMATIC ERRORS 



In this section we estimate the systematic uncertainties in the S'f/(3)-breaking ratios 
IbsI Ib^ ^- For clarity, we present each source of error in a separate subsection. The 
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total error budgets for both quantities are given at the end of the section in Table I Villi 



A. Chiral extrapolation fit ansatz 

As described in the previous section, we extrapolate our data to the physical quark masses 
using a linear fit matched onto NLO SU(2) HM^PT at the value of our lightest data point. 
The SU (2) expressions are derived using the symmetries of the lattice theory; therefore they 
contain the correct dependence of the i?-meson decay constants and mixing matrix elements 
on the quark mass and lattice spacing through NLO when the quark masses are sufficiently 
light. Given the large light-quark masses used in our simulations, however, it may be that 
the SU{2) chiral logarithms do not become important until the pion mass is even lighter 
than the range of our data. Thus we vary the location of the matching point as one way 
to estimate the systematic uncertainty due to the chiral extrapolation. In addition, because 
we use NLO SU{2) HMxPT to obtain our central value, we must estimate the systematic 
uncertainty due to the truncation of higher-orders in HM^FT. We do this in several ways: 
(i) by explicitly adding higher-order terms to the linear plus SU{2) HMxPT fit function, (ii) 
by matching the linear fit onto the NLO SU{3) HMxPT expressions, and (iii) by varying 
the value of the low-energy constant / in the coefficient of the NLO chiral logarithms. 

Because we do not know a priori at what mass the SU(2) chiral logarithms become 
important, we vary the point at which we match the linear fit onto NLO SU{2) HM^PT 
in order to estimate the systematic uncertainty due to the choice of matching point. At 
the matching point used in the preferred fit, the ratio of the light up-down quark mass to 
the strange valence-quark mass is ami/amg ~ 0.28 and to the strange sea quark mass is 
ami/amh = 0.25. Since both of these quantities are small, we expect the strange quark can 
be integrated out of the chiral effective theory and that SU (2) HMxPT is applicable in this 
region. As the light-quark mass decreases, SU(2) becomes an even better approximation. 
At light-quark masses above this point, however, the ratios ami/ arris and ami/amh are no 
longer small expansion parameters, and SU{2) will eventually cease to apply. Therefore, 
when we estimate the systematic error due to the choice of matching point, we only consider 
extrapolations in which the matching point is closer to the chiral limit than in the preferred 
fit. The limiting case is where the matching point is at the chiral limit, which corresponds 
to a purely linear extrapolation. Although the choice of a linear fit is not based on effective 

31 



field theory, it cannot be ruled out by the data. We therefore take the difference between 
the extrapolated values for fBs/fsa ^^'^ ^ obtained with the preferred fit and with a purely 
linear fit to estimate the systematic error due to the choice of matching point; this leads to 
an error in fsj fsa of 6.6% for both the APE and HYP data and in ^ of 6.8% (6.9%) for the 
APE (HYP) data. We compare the results of the preferred fit with those of the purely linear 
extrapolation in Fig. [7] for ^Bs/^Bi and in Fig. [8] for Comparison with the 

linear extrapolation leads to a conservative systematic error estimate since we know that, 
when the pion mass is sufficiently light, the chiral logarithms will become important and the 
extrapolation function will begin to curve upward. However, the linear fit provides a clear 
lower-bound on the possible extrapolated value. 



_ 'linear + SU[2) HMxPT: *ij,/*B,, = 1.165(i 
if^lme-Ai: *B,/*fl^ = 1.088(82) 
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linear + SU{2) HMxPT; *b,/*b^ = 1. 153(47) 
jg^lmear: $b^/$b, = 1.077(44) 
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FIG. 7: Chiral extrapolation of '^Bs/'^Bi = \/'^BsI'^Bi ■ Ibs/Ibi using a linear fit matched onto 
NLO SU{2) HMxPT at the lightest data point (solid band) and a linear extrapolation all the 
way to the chiral limit (hatched band). Errors shown are statistical only. The left plot shows 
the comparison for the APE data, while the right plot shows the same comparison for they HYP- 
smeared data. 



In order to estimate the error due to the omission of higher-order terms in the chiral 
expansion, we compare the result of the linear plus NLO SU (2) HM^PT fit to a fit sup- 
plemented by NNLO analytic terms. Once the valence and sea strange quarks have been 
integrated out of the chiral effective theory, the SU{2) HMxPT expressions can only depend 
upon the light up-down quark mass mi and the lattice spacing a. Thus the only possible 
NNLO analytic terms are those proportional to mf, mia"^, and a^. Because we are working at 
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FIG. 8: Chiral extrapolation of a/ MbJMbi = \frriBjrriB^ {fss \/Bbs/ fsi \/Bbi) using a linear fit 
matched onto NLO SU (2) HM^PT at the lightest data point (solid band) and a linear extrapolation 
all the way to the chiral limit (hatched band) . Errors shown are statistical only. The left plot shows 
the comparison for the APE data, while the right plot shows the same comparison for they HYP- 
smeared data. 

a single lattice spacing, the contribution proportional to is indistinguishable from the LO 
normalization factors R<s, and Rm in Eqs. ( ITHj) and (179|) and the contribution proportional 
to m/a^ is indistinguishable from the NLO terms Cimi and Dimi. Thus there is only one ad- 
ditional free parameter at NNLO. In practice, we implement the "NNLO fit" by performing 
a quadratic fit in mi to our data and then matching to the NLO SU (2) HMxPT expressions 
supplemented by the NNLO analytic term proportional to mf at the value of our lightest 
data point. When we include the term quadratic in mi, however, the extrapolated values 
for fBs/fBd ^ have significantly larger statistical errors (about 20 — 40%) than those 
obtained with the preferred fit (about 4 — 8%). Therefore we cannot draw any meaningful 
conclusion about the change in the central values of fss/fBa and ^, since the outcome of the 
NNLO extrapolation is consistent with the NLO extrapolation within the large statistical 
errors. This increase in statistical errors is to be expected since we have introduced an 
extra free parameter, and it is difficult to constrain two parameters with only three data 
points. Thus we do not use the "NNLO fit" as a way to estimate the chiral extrapolation 
error in this work, but leave it as a way to estimate the uncertainty due to the omission of 
higher-order terms in future analyses when we have more data points and smaller statistical 
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errors. 

We can also estimate the systematic uncertainty in the chiral extrapolation by comparing 
the values of /b^//_b^ and ^ obtained by matching onto SU{2) HM^PT with those obtained 
by matching onto SU{3) HM^PT. Although the SU{2) chiral logarithms are a subset of the 
SU (3) chiral logarithms, the SU (2) and SU (3) theories have different series expansions and 
different degrees-of-convergence within their ranges of applicability. Therefore the compar- 
ison with NLO SU{3) HM^PT provides another means of estimating the error due to the 
use of NLO SU{2) HM^PT. The expressions for ^[/(3)-breaking ratios at NLO in SU{3) 
HM^PT are, schematically, 

$B 2B _ 

" - 1 + "chiral logs" + „.„ Cvai(m^ - mi), (80) 



^B, " (4vr/) 



^ = 1 + "chiral logs" + j^^d^m, - mi), (81) 

where "chiral logs" indicate non-analytic functions of the pseudo-Goldstone meson masses, 
e.g. m| log(m|/A^). These are derived from Eqs. ( 1D2|) and ( 1D3[) by taking the ratio of the 
expressions for the valence quark x = s over the expressions for x = I. Because the strange 
quark is treated in the same manner as the up and down quarks in the SU (3) chiral effective 
theory, the low-energy constants are the same in the numerator and denominator. Thus 
the overall normalizations cancel in the expressions in Eqs. ( IHOj) and (IHT]) and the ratios 
are constrained to unity in the limit rrii — j- m^. Therefore the expressions for the SU{3)- 
breaking ratios have one free parameter each, instead of two as in the SU{2) HM^PT case. 
This means that the SU (3) HM^PT expressions cannot be matched smoothly onto the linear 
fit of the heavy data. We choose to make the value of the extrapolation function continuous, 
while leaving a discontinuity in the slope at the matching point. The difference between the 
linear plus NLO SU{2) HM^PT fit and the linear plus NLO SU{3) HM^PT fit leads to a 
difference in JeJIb, of 2.3% (2.4%) for the APE (HYP) data and in ^ of 2.6% (2.5%) for 
APE (HYP). 

In the preferred linear plus NLO SU (2) HM^PT fit, we set the leading-order pseudoscalar 
meson decay constant / equal to the experimentally- measured value of This fixes the 



coefficient of the chiral logarithms and improves the convergence of the chiral expansion 
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64( 1 . At NLO in xPT, however, it is equally consistent to use the pseudoscalar decay constant 



in the SU{2) chiral limit /o or the kaon decay constant fx because the different choices only 
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affect NNLO terms tliat are of fiigher-order tlian we consider. We tlierefore vary / between 



/o = 115 MeV 



34| and fx = 155.5 MeV {g] in order to estimate tlie systematic uncertainty 
due to the omission of liigher-order terms in the chiral expansion. The use of /o in the SU{2) 
chiral hmit, which the RBC and UKQCD Collaborations find to be about 10% lower than 
/tt, leads to a difference from the central value for fsJlBd of 2.1% (2.2%) for APE (HYP) 
and from the central value for ^ of 2.3% for both APE and HYP. The use of Jk leads to 
similar changes in fss/ fsa of 2.1% and in of 2.2%. 

We take the largest of the uncertainties enumerated above, which is obtained from the 
difference between the preferred linear plus NLO SU{2) HM^PT fit and the purely linear 
fit, for the final estimate of the chiral extrapolation error; this leads to the values in the row 
labeled "chiral extrapolation" in Tab. IVHIl 

B. Uncertainty due to gB*BiT 

Although we fix the value of the B*-B-'k coupling (and hence the coefficient of the one- 
loop chiral logarithms) in the chiral extrapolation of our lattice data, gB*Bn is in fact poorly 
known from phenomenology. Hence we must vary the value of Qb^b-k over a sensible range 
based on lattice QCD calculations and phenomenology in order to estimate the systematic 
errors in fss/fBj, and ^ due to the uncertainty in the coefficient of the one- loop chiral 
logarithms. 

There has only been one unquenched lattice QCD calculation of the B*-B-tx coupling 
with a complete associated error budget, which give s qB *BTv = 0.516(5)(50), where the first 



uncertainty is statistical and the second is systematic 62|]. Although this value was computed 
in the static heavy-quark limit and neglects the effects of the dynamical strange quark, we 
take this as our central value because of its small statistical errors and full systematic error 
budget, as well as because it lies in the middle of the range of values presented in the 
literature. Another recent determination of Qb'Btt = 0.44 ± 0.03tooo in the static limit in 
two-flavor lattice QCD, where the errors are due to statistics and the chiral extrapolation 



uncertainty, is consistent with this result 



65|. QCD sum rules and the relativistic quark 

The value of the B*-B-7r 



model give a slightly lower value of gB*Bn ~ 0.38 ± 0.08 
coupling is expected to be close to the D*-D-ti coupling because of heavy-quark symmetry. 
We can therefore also use lattice QCD calculations and phenomenological extractions of 
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Qd* Dn as estimates. Becirevic and Haas recently computed gD*Dn = 0.71(7) in 2-flavor lattice 
QCD, but this result is from only a single lattice spacing and presents no estimate of the 



data 
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systematic error |67|. Stewart computed the value of gD*D-K in 1998 by fitting to experimental 



and recently updated his result to include the experimental measurement of the 
D* decay width |^|. His latest determination is Qo'D-k = 0.51, but with no error quoted 70 |. 



The most sophisticated extraction of the D*-D-7t coupling was performed by Kamenik and 
Fajfer, and gives go'Dn = 0.66lQ;Qg, where the uncertainty only reflects the error due to 
counterterms 71|. Finally, we note that the chiral-continuum extrapolations of 5-meson 
quantities by the Fermilab/MILC Collaboration 72| and the HPQCD Collaboration |l^ 



tend to prefer even smaller values of the B*-B-tx coupling than those in the literature. 

Most of the results presented above do not have complete error budgets, and are incon- 
sistent within the quoted errors, so for this work we take Qe'E-k = 0.516 ± 0.2 to account 
for the spread of values. We then vary gB*B-K within this range to determine how much it 
changes the central values for fssl Ib^ and ^. This leads to an uncertainty in fsJIsa of 
3.2% for both APE and HYP and in ^ of 2.1%. 

C. Discretization errors 

In this work we only analyze data at a single, relatively coarse lattice spacing of 
a ~ 0.11 fm, so we estimate the size of discretization errors with power-counting. As a 
consistency-check of our estimation procedure, however, we can compare the estimated er- 
rors in the individual decay constants and matrix elements with the observed differences 
in those quantities for the APE and HYP-smeared data. This is because, aside from sta- 
tistical errors, the differences in the values obtained from the two smearings are due to 
discretization effects and higher-order corrections to the renormalization factors. We find 
about 15-20% differences in the decay constants and matrix elements obtained from the 
APE and HYP-smeared data, whereas we estimate by power-counting that discretization 
errors in fs^ should be about 15% and in JUg should be about 20%. Thus our observations 
are consistent with the scaling behavior expected from power-counting. We can also compare 
our power-counting estimates with the findings of the ALPHA Collaboration, who performed 
a study of Jb, at several lattice spacings using static 6-quarks with similar link-smearings 



in the quenched approximation 73|. Although ALPHA observes a violation of O(a^) scal- 
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ing behavior at inverse lattice spacings below ~ 2.5 GeV, the difference between the 
predicted value for /b, at a ^ 0.11 fm given C(a^) scaling and the value of fs^ that they 
obtain in the continuum limit using data within the scaling region is only about 25%, which 
is again close to our power-counting estimate. Thus we expect that naive power-counting 
should lead to a reasonable estimate of the discretization error in the SU (3)-breaking ratios. 

For the error estimates in this subsection, we evaluate the strong coupling constant at 
the lattice scale, af^{l/a) ~ 1/3. We choose Aqcd = 500 MeV because the typical QCD 
scale that enters heavy-light quantities tends to be larger than for light-light quantities, as 



indicated by fits to moments of inclusive 5-decays using the heavy-quark expansion [74 1. 
Fortunately, some of the finite lattice-spacing effects cancel in the ratios fss/fBd ^^^^ ^- This 
can be seen by the fact that, although the APE and HYP data differ by about 15-20% for 
the individual decay constants and matrix elements, they agree within statistical errors for 
the ratios. In S'f/(3)-breaking quantities, errors must be proportional to the difference in 
quark masses {rris — m^. Dimensional analysis therefore suggests that contributions to the 
total discretization error are suppressed by the factor {fhs — ^d)/AQCD ~ 1/5, where we 
use nis and rhd to denote the renormalized quark masses in the MS scheme p] (as opposed 
to the bare lattice quark masses) in this subsection and the next. The observed size of 
S'f/(3)-breaking effects in the i?-meson decay constants (/b^//b^ — 1) and in the 5— mixing 
matrix elements (^ — 1) are consistent with this expectation. Discretization errors in fss/fB^ 
and ^ can arise from both the actions and the operators. We estimate each source of error 
separately, and add them in quadrature to obtain the total discretization error. 

None of the actions that we are using are (9(a^)-improved. Therefore the leading dis- 
cretization errors from the domain-wall fermion action and Iwasaki gauge action are of 
(9(a^AQQj-,). When combined with the 5'f/(3)-breaking suppression factor, this leads to dis- 
cretization errors in the ratios fss/ ^ of ^('^^^qcd ^ — f^d) / Aqcd) ~ 1-7%. The 
leading heavy-quark discretization errors from the static action are also of Ola'^A^Qj^). Hence 
heavy-quark discretization errors also contribute ~ 1.7% to the total error in the ratios. Be- 
cause we improve the heavy-light axial current used to compute the decay constant through 
0{asap), the leading discretization errors from the heavy-light current are of 0{asam), 
0{a'j.a Aqcd), and 0{a'^AQQj)).^ When combined with the S'f/(3)-breaking suppression fac- 



There are also discretization errors from mixing with operators of other chirahties that are proportional 
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tor, this leads to discretization errors in the ratio fssl fs^ of 0{as x (am^ — arrid)) ~ 1.2% 
plus 0{alaAQCD x (m^ - md)/AQCD) ~ 0.6% plus (^(a^A^cD x (m^ - md)/AQCD) ~ 1.7%. 
Although we do not improve the heavy-light four fermion operator used to compute the 
fi-mixing matrix element, the operator does not have any tree-level 0{a) errors 



Thus the leading discretization errors in the ratio ^ from the four-fermion operator are 
of 0{as X {anis - anid)) ~ 1.2% plus O^asaKqcv) x [nis - md)/AQCD) ~ 1.9% plus 
C(a2A|cD X (^s - md)/AQCD) ~ 1.7%. 

Adding the contributions from light-quark and gluon discretization errors, heavy-quark 
discretization errors, and discretization errors in the heavy-light current or four-fermion 
operator in quadrature, we estimate the error in fsJfBa to be ~ 3.2% and the error in ^ to 
be ~ 3.7%. 



D. Heavy-light current and four-fermion operator renormalization 

We compute the renormalization factors needed to match the lattice axial current and 
four-fermion operator to the continuum using one-loop lattice perturbation theory. This 
leaves a residual error due to the omission of higher-order terms. Based on power-counting, 
we estimate the truncation error in the coefficients to be of 0{al), which is the size of the 
first neglected term in the series. As we noted earlier in Sec. IIII C[ however, the matching 
coefficient Z$ cancels in the ratio of decay constants fss/ fsai thus its contribution to the 
error in fss/fBa is zero. Although such an exact cancellation does not occur for the ratio 
of mixing matrix elements ^, the error in due to the uncertainty in the ratio of matching 
coefficients Zsp/Zva is suppressed by the S't/(3)-breaking factor {fhs — ma)/ Aqcb- This is 
because, in the SU{3) limit, the four-fermion operator matrix elements would be equal in the 
numerator and denominator, so the error in ^ from the renormalization factor uncertainty 
would be zero. We therefore expect the error in fssl fs^ to be 0% and the error in ^ to be 
of 0{al X {fhs — w^d)/AQCD) ~ 2.2%. This error will decrease with the inclusion of data at 
a finer lattice spacing because the smaller coupling constant will improve the convergence 
of the series. 



to amies- These effects, however, are expected to be sub-percent level in the matrix elements 
therefore negligible in the S'J7(3)-breaking ratios. 



751, and 
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E. Scale uncertainty 



Because in this paper we only compute dimensionless ratios, the uncertainty in the de- 
termination of the lattice spacing only enters implicitly through the uncertainty in the 
light-quark masses and in the renormalization factors. We estimate the systematic error due 
to the truncation of lattice perturbation theory in the previous subsection and due to the 
light-quark mass determinations in the following subsection. 



F. Light- and strange-quark mass uncertainties 

We obtain the physical decay constants and mixing matrix elements by setting the light- 
quark masses to their physical values in the linear plus SU (2) HMxPT chiral extrapolation 
formulae, once the low-energy constants have been determined from fits to numerical lattice 
data. We use the bare-quark mass value determined from fits to the light pseudoscalar 
meson masses 34l |: 



am^rf + am^es = 0.001300(85), (82) 

where rriud is the average of the up and down quark masses. The quoted error includes 
both statistics and the systematic uncertainties from the chiral extrapolation, finite-volume 
effects, discretization effects, and the unphysical strange sea quark mass. In order to estimate 
the systematic uncertainty in the ratios fss/fBa and we vary the bare light-quark masses 
within their stated uncertainties. We then take the maximal difference from the central value 
to be the systematic error. From this method, we find that the systematic error in fsJlBa 
due to the uncertainty in the light-quark mass determination is 0.2% for both smearings and 
the systematic error in ^ is 0.2% regardless of the smearing used for the heavy quark. 

Because the strange-quark mass does not explicitly appear in the SU (2) HM^PT extrap- 
olation formulae, we cannot estimate the errors in fss/fBa and ^ due to the simulated valence 
and sea strange-quark masses with the simple method used above for the light quarks. We 
must instead address the errors due to the uncertainty in the strange valence-quark mass 
(which is set to the physical value of m^) separately from those due to the choice of strange 
sea-quark mass (which is not tuned to the physical m^). 

We calculate the decay constants and mixing matrix elements with the strange valence- 
quark mass tuned to the physical value obtained from fits of the light pseudoscalar meson 
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masses 



19| 



am, + arrives = 0.0390(21), (83) 

where the error includes those due to statistics and to uncertainties in the pseudoscalar meson 
masses, residual mass, and lattice scale. Since the time at which our numerical computation 
was performed, however, we have improved the strange-quark mass determination using 



data on the larger-volume 24^ ensemble. This analysis yields arus + amj.es = 0.0375(17) 34j . 
where the error includes statistics and the systematic uncertainties from the chiral extrap- 
olation, finite-volume effects, discretization effects, and the unphysical strange sea quark 
mass. Fortunately, the two determinations are consistent within their stated uncertainties 
and the small change in valence strange-quark mass leads to a negligible difference in the 
the decay constants and mixing matrix elements. A linear interpolation of and Mb^ 
from arUs = 0.0359 0.0343 leads to at most a percent-level change in these quantities, 
which is too small to resolve within our large statistical errors. We therefore conclude that 
it is sufficient to forgo the interpolation to arris = 0.0343 in our current analysis. 

Although the strange-quark mass is integrated out in SU(2) HMxPT, the value of the 
strange valence-quark mass enters the chiral extrapolation described in Sec. |V]via the linear 
fit in the heavy mass region that is constrained to unity in the limit mi — j- m^. Once 
the slope in the linear region has been determined from fits to numerical lattice data, we 
estimate the errors in and ^ due to the uncertainty in by fixing the slope to 

this value and varying m^ within its stated uncertainty given in Eq. (!83|) . This shifts 
the location of the SU{3) limit and hence the linear portion of the extrapolation in the 
horizontal direction by an amount ±crm^; consequently it shifts the location of the point 
at which we match onto SU{2) HMxPT in the vertical direction by the slope times ±0"^^. 
From this method, we find that the systematic error in fsJfBa due to the uncertainty in 
the valence strange-quark mass determination is 0.4% regardless of the smearing used for 
the heavy quark and the systematic error in C, is 0.4% (0.3%) for the APE (HYP) data. 

The value of the strange sea-quark mass does not explicitly enter the chiral extrapolation 
formulae. Nevertheless, it implicitly affects the determinations of and ^ through 

the values of the SU(2) HMyPT low-energy constants. The latest results of the RBC and 

n n 

UKQCD Collaborations for /^r, fx, and Bk [76| rely on reweighting [77| to self-consistently 
tune the strange sea-quark mass to its physical value during the chiral extrapolation of the 
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pseudoscalar meson masses and decay constants. We can therefore estimate the effect of 
the incorrect strange sea-quark mass on fssl fsa and ^ based on observations in the hght 
pseudoscalar meson sector. Prehminary studies by RBC and UKQCD find that reweighting 
to the physical strange sea-quark mass leads to no statistically significant change in 771^^ 



and a 2% change in /^r [78|. If we assume that the unphysical strange sea-quark mass 
leads to a 2% error in the decay constants and matrix elements, then it will lead to a 
2% X {rhs — w^d)/AQCD ~ 0.4% error in the 5'f/(3)-breaking ratios fss/ Ib^ and C,. 

In order to obtain the total systematic error in fssl fsd and ^ due to the light- and 
strange-quark masses, we add the three estimates for the light-quark error, the strange 
valence-quark error, and the strange sea-quark error in quadrature. This leads to an error 
in fsJlB^ of 0.6% for both the APE and HYP data and an error in ^ of 0.6% (0.5%) for 
the APE (HYP) data. 



G. Finite-volume errors 



We estimate the uncertainties in the ratios /^^ / /^^ and ^ due to finite volume effects using 
one-loop finite-volume HM^PT. The effect of the finite spatial lattice volume is simply to 
turn the one-loop integrals that appear in the HM^PT expressions into sums, which we 
show in Eqs. flD13p - (lD14p . It is therefore straightforward to compute the corrections to 
our data given the parameters of our simulations. We find that the one-loop finite volume 
corrections are about 1% at our lightest quark mass and even less at heavier masses. We 
therefore take 1% to be the systematic uncertainty in fss/fBa and due to the finite spatial 
lattice volume. 



H. l/rrif, corrections 

Because we work in the static heavy-quark limit, our results for JbsI fsd and neglect 
relativistic effects due to the finite fe-quark mass. We therefore estimate the size of the 
omitted relativistic corrections using power-counting. The leading l/nib corrections to the 
decay constants and matrix elements are of 0{KQCY)/m}j). Because we are only computing 
the S't/(3)-breaking ratios fsJfB^ and ^, however, relativistic effects in these quantities are 
further suppressed by a factor of (m^ — rnd) / Aqcd- Therefore we expect the relativistic 
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TABLE VIII: Total error budget for the 5C/(3)-breaking ratios fss/fBa Each source of 

uncertainty is discussed in Sec. IVIl and is rounded to the nearest percentage. 



fBs/fsa ^ 



uncertainty 






A V>'f? 


rl 1 r 


statistics 


o/o 


A 07 


aW 
D70 


070 


chiral extrapolation 


7% 


7% 


7% 


7% 


uncertainty in gs* Bw 


3% 


3% 


2% 


2% 


discretization error 


3% 


3% 


4% 


4% 


renormalization factors 


0% 


0% 


2% 


2% 


scale and quark mass uncertainties 


1% 


1% 


1% 


1% 


finite volume error 


1% 


1% 


1% 


1% 


l/mb corrections 


2% 


2% 


2% 


2% 


total systematics 


9% 


9% 


9% 


9% 



corrections to these quantities to be approximately 



^QCD ^ 'Tl, 



A 



2%, 



QCD 



(84) 



where the scale Aqcd cancels in the ratio, and we use the average MS quark masses listed 



in the PDG 



6|. 



VII. RESULTS AND CONCLUSIONS 



Using the experimentally-measured ratio of masses m^o/m^o = 5366.6/5279.5 



1.0165 



6[ , we obtain the following values for the SU (3)-breaking ratios of 5- meson decay 



constants and mixing matrix elements: 

f 1.16(09)(10) APE f 1.13(07)(10) APE 

— ^ = <J and ^ = <J , (85) 

[ 1.14(05)(10) HYP [ 1.13(05)(10) HYP 

where the first errors are statistical and the second are the sum of the individual systematic 
errors added in quadrature. We find good agreement between the different link-smearings, 
indicating that, despite the use of a single lattice spacing, discretization errors are small 
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in the ratios. We therefore average the APE and HYP determinations to obtain our final 
results. After adding the statistical and systematic errors for each link smearing in quadra- 
ture, we compute the average assuming that the two determinations are 100% correlated 



using the method of Ref. 



79j: 



^ = 1.15(12) (86) 



e = 1.13(12), (87) 

where the errors refiect the combined statistical and systematic uncertainties. Although we 
computed these quantities in the static 6-quark limit, the inclusion of the neglected l/rrib 
corrections (which we estimate to be about 2%) produces a negligible change in the total 
errors in Eqs. (15^ and fl87p given the size of our other uncertaintites; thus our results can be 
directly compared to phenomenological determinations and other lattice QCD results using 
relativistic 6-quarks. As shown in Fig. [9l our results agree with the published results of the 
HPQCD Collaboration (^ = 1.258 ± 0.025stat. ± 0.021sys,) Q and the preliminary results of 



the Fermilab Lattice and MILC Collaborations (^ = 1.205 ± 0.037stat. ± 0.034sys) [15|. 

Although our results have significantly larger errors than the other Nf = 2 + 1 fiavor 
determinations, in this work we have demonstrated the viability of our lattice computation 
method. In particular, we have introduced the new approach of using SU{2) heavy-light 
meson chiral perturbation theory to extrapolate Nf = 2 + 1 lattice QCD results for 5-meson 
quantities to the physical quark masses. The largest sources of error in our calculation are 
from statistics and the chiral extrapolation, and we expect to reduce the sizes of both in a 



future work that analyzes the 24^ domain- wall ensembles with the same lattice spacing 34 1. 
Some of the 24^ ensembles contain almost three times as many configurations as we have 
analyzed in this work. Furthermore, the use of a larger spatial volume will allow us to 
simulate at lighter valence and sea quark masses. Once we have made these improvements, 
our results will provide a valuable cross-check of these important inputs to the CKM unitarity 
triangle analysis and determination of the ratio of CKM matrix elements |Vtd|/|Vts|. 
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Appendix A: SU{3) projection 



In this appendix we describe the two SU{3) projection schemes used in this paper and 
show their equivalence in the weak couphng hmit. Necessary properties of the projector 
mapping an arbitrary, complex 3x3 matrix on the SU{3) subgroup are idempotence and 
gauge-covariance. These properties are not sufficient to specify a unique projector and 
hence several choices exist. For APE smearing we use the unit circle projection which 
is based on polar decomposition, while for HYP smearing we seek iteratively the matrix 
f^max £ SU{3) which maximizes Re Tr(f/jnax^^) 0]- 



First we describe the unit circle projection [46[. For a complex 3x3 matrix V with 
det(V^) 7^ 0, we calculate the matrix 

W = V[V^V]-'^\ (Al) 

which is unitary by construction and has a spectrum lying on the unit circle. The square 
root is obtained by Jacobi matrix diagonalization. From W we obtain a special unitary 
matrix by computing 

V = [{det{W)]-^/^W. (A2) 

This projection is idempotent since an element of SU{3) is projected by Eqs. (1A1|) and f lA2p 
back onto itself. The projection is also gauge covariant, as we now show. The matrices V 
and V^^ transform as 

V Gl^G^j and GrV^G[, (A3) 

and hence 

V^V GjjV^VgJj. (A4) 

Since [V^V]-'^/^ has the same transformation properties as V'^V one finds for the transfor- 
mation of W 

W GiWC^ji, (A5) 
from which the gauge-covariance of V follows. 
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Next we describe the projection method used for HYP smearing. This method is based 



on seeking a matrix f/max for which 



Umax G S'[/(3) iRe Tr([/max^^) IS maximal . 



(A6) 



The maximizing matrix^ [/max 



as outhned in Refs. 
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is found iteratively by decomposing it into its SU (2) subgroups 
80|. Although, by the compactness of SU{3), a global maximum 



exists, it is not guaranteed that this iterative procedure converges to the true maximum. 
Furthermore, the matrix obtained by the iterative procedure depends on the details of the 
iteration algorithm. In practice, however, one finds that f/max is numerically very close to 
V, the SU{3) projection of the matrix V obtained from the unit circle projection, 46|, and 
hence the iteratively found maximum is close to the true maximum. Thus we choose 



V = U„ 



(A7) 



as the SU{3) projection of the complex matrix V . This projection is idempotent because for 
V G SU{3), unitarity and the triangle inequality imply the bound Re Tr(yi^"'') < 3 which 
is uniquely saturated hj V = V. By construction this projection is also gauge-covariant. 

Finally, we show that although the two projection methods are generally not equivalent 
they differ only at second order. For this discussion we assume that the iterative method 
converges to the global maximum and that the unprojected matrix V is approximately 
unitary and unimodular (which is true in the weak coupling limit). 

In order to compare the two projection schemes, we first write the polar decomposition 
of V as 

V = e'^VH, (A8) 

where 6 G (— 7r/3, 7r/3], if is a Hermitian matrix with eigenvalues Ai,A2,A3 > 0, and V is 
the S'f/(3)-projected matrix given by Eq. flA2p . We then diagonalize H = MAM'^ [A = 
diag(Ai, A2, A3) and M G SU{3)] and substitute Eq. (jAll) into Eq. (jAe])- The resulting 
function to be maximized is 



ReTr(e-^^f/'A) = Re(e-^''(Ain;i + Aa^s + A3n'33)), 



(A9) 



where u[j is an element of U' = M^V^UmaxM G SU{3). If the maximization is carried out 
over f/(3) rather than SU{3), we could choose u[-^ = u'22 = M33 = £ 



This satisfies the 
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bound Re Tr(e *^[/'A) derived from the triangle inequality. If we maximize over SU{3) and 
V is already an element of SU{3) then it follows U' = 1 and there exists a unique maximum 
for 6* = such that V = V. 

For the generic case maximizing over SU{3), the constraint det(t/') = 1 forces a maxi- 
mum at f/' 7^ 1, and hence V ^ V. Nevertheless, when V is approximately unitary and 
unimodular we can write 6 = 66 with \66\ ^ 1. Then Afc = 1 + 6Xk and \6Xk\ <^ 1. From 
that follows u'),). = 1 + i6k with ^ 1 {k = 1, 2, 3), and we can work to first order in all 
"small" quantities. Although U' can in principle develop small, off-diagonal elements, the 
constraint U' G SU (3) always imphes that 5/; G M and 6i + 62 + = 0. Therefore we must 
maximize 

3 

^ Re (1 + 6Xk + t6k - 166) , (AlO) 

k=l 

subject to the constraint 81 + 82 + 5^ = 0. Because this is independent of 66, this is equivalent 
to the case 6 = discussed above. Consequently, to first order V = V. This explains why 



the two projection schemes lead to numerically close results as reported in 



46|. 



Appendix B: Continuum QCD to HQET matching and HQET running 

In Eqs. flMl) and fl32p we separated the matching coefficients of the operators into three 
contributions: the QCD to HQET matching factor at scale ruh, the HQET running from rrib 
to rric, and the HQET running from to fi. Here we collect the results for these factors 
and present them for general numbers of flavors Nf and colors N^. The Casimir factors that 
appear in the expressions are 

Cf = ^^, Ca = N,, and = i. (Bl) 

The one-loop QCD to HQET matching factor for the axial current operator renormalized 
in the MS scheme with naive dimensional regularization is calculated in Refs. 
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C^K) = 1 - 2C^^^ + 0{al), (B2) 



while the renormalization group (RG) evolution factors entering Eq. (13T|) are given by 48l |49 1 
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The coefficients in the evolution factor are 



d' 



J 



{Nf) 



7?V2/3o, 



(B4) 
(B5) 



with 



/3o 
/Si 



3 

34 



4 



20 



CI - iCpTpNf - —CaTfNj, 



3 ^ ... 3 
and the one- and two-loop anomalous dimensions 



(B6) 
(B7) 



if 
7i" 



-16Cf 



49 
96 



Ca 



32 



24' 



(B8) 
(B9) 



The one-loop QCD to HQET matching factors for the AB = 2 four quark operator are 



calculated in Refs. 



50 



53|: 



Zi(mb,mb) 



8N^ + 9Nc-15 as{mb) 



2Nr 



An 



An 



(BIO) 
(Bll) 



The operators 0^ and Os mix under RG evolution (see Eq. (!32l) ). such that the evolution 
factors for the matching coefficients can be written as a 2 x 2 matrix: 



U 



Note that the matrix element 1/^2 ^ '^ is zero to all-orders in perturbation theory due to the 
fact that heavy-quark spin symmetry prohibits mixing from Os into O^. Because Z2 in 
Eq. fIBlip has no tree- level contribution, the two-loop expression for f/|f^^ and one- loop 
expressions for ^72^^'' ^22^'^ sufficient to determine the 0{as) matching coefficients. 



"^11 







(B12) 



22 
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The matrix elements in Eq. f lB12|) are given in Refs. 5l|-l53|: 



with the coefficients 



V 



hOK^), (B13) 
+ 0{as), (B14) 

(B15) 



and the anomalous dimensions 



(0) 



7 



1 + 1/iV, 



-6Cf + 4('1 



iVc-1 

12Nr 



(^127A^,2 + 143iVc + 63 - 57 /N^ 



+87r\N^ - 2N^ + A/N,) - Nf{28Nc + 44)) . 



(B16) 
(BIT) 



(B18) 



(B19) 



Appendix C: Finding matrix elements from Green's functions 

In this appendix we provide a derivation of Eqs. f l63p and fl72p which give the relation- 
ship between the Euclidean space correlation functions which are the direct results of our 
simulations and the matrix elements between normalized Hilbert space states which we wish 
to determine. In this appendix, all quantities are written in lattice units. The derivation 
of Eq. f l63|) follows the standard steps of inserting a complete set of energy eigenstates into 
the Green's functions C^^(t,t()) and C^^(t,to) which are defined in Eqs. ( j59l) and (j60l) and 
appear in the numerator and denominator of Eq. ( l63l) : 

C^^(t,to) = $^$^(0|Ao"(f,0)|n)(nK(0)t|0)e-^"(*-*°) (CI) 

n x€V 

C^^{t,to) = ^|(n|AS^(0)t|0)|V^"(*-*o\ (C2) 

n 
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where £"„ is the energy of the state |n), a state with unit normahzation. Because of the 
translational invariance of either the wall source operator or the sum of the local operator 
over a temporal hyperplane, the intermediate state |n) must have zero momentum. There 
is then a single, finite- volume state, \Bq{p = 0)v), separated from the other excited states 
by a non-zero energy gap, which will dominate as the time separation t — to becomes large. 
Thus, in the limit of large t — to, Eqs. (ICip and (1020 become: 

C^^(t,to) = L^O\A^,mB,{p = 0)v){B,{p = 0)vK(0)t|0)e-"^^.(*-*«) (C3) 
C^^(t,to) = |(i?,(p = 0)yK(0)t|0)|2e-"*^.(*-*°\ (C4) 

where we have used translational symmetry to replace the sum over x in Eq. (ICip with the 
factor L^. 

Finally, Eq. f l63p follows by taking the ratio of Eq. (1C3P and the square root of Eq. (IC4p 
to remove the unwanted factor (OI^q (0)|-Bq(p = 0)) and recognizing that in the limit of large 
volume the unit-normalized state \Bq{p = 0)y) in Eq. fICSP and the covariant, delta-function 
normalized state \Bq{p = 0)) in Eq. (!57|) are related by: 

^2mBL^\Bq{p = 0)y) ^ \Bq{p = 0)). (C5) 

This discussion is standard and has been repeated here to provide a familiar background 
for the derivation of Eq. (1721) where a new approach, special to the static approximation, 
is required. The complication in Eq. fl72|) arises because of the use of a box source which 
is not translationally invariant and which produces a superposition of i?-meson states with 
various momenta. In the static limit, all of these states are degenerate since their energy 
no longer depends on their momenta. Thus, we cannot assume that the large time limit, 
t f ^ t ^ to will project onto a unique ground state. However, as is worked out below, we 
can use an additional symmetry of the static approximation to show that the normalization 
of the box source cancels between the numerator and denominator of Eq. (172|) js?! . 

Fortunately, this large set of degenerate states resulting from the momentum indepen- 
dence of the energy of the heavy-light meson, can be distinguished by a new conservation 
law which becomes exact in the static limit: the local conservation of heavy-quark num- 
ber. The absence of spatial derivatives in heavy-quark actions shown in Eqs. ( lT2l) and (lT6ll 
implies that the total number of heavy quarks at each spatial site is separately conserved. 
This conservation law results from the invariance of the heavy-quark action under the phase 
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rotation of the heavy-quark Grassmann variables: 

h{x) e'^^'^^hix) (C6) 
h{x) e~'^^^^h{x) (C7) 

in which a different phase 6{x) can be used for each spatial site x. In the quantum mechanical 
Hilbert space this symmetry corresponds to a family of operators Nh{x) = h{x)'~f^h{x) which 
commute with the Hamiltonian (continuum theory) or the transfer matrix (lattice theory). 
Thus, when interpreting our Green's functions we can introduce a complete set of S-meson 
states \Bq{x)) which are eigenstates of both the Hamiltonian and the number operators 
Nh{x') with a standard normalization: 

Nh{x')\B,{x)) = -6^,x'\B,ix)) (C8) 
{Bg{x)\B,{x')) = 6,,,,. (C9) 

Here the Kronecker delta, 5^,^', used above is appropriate for the discrete positions that 
appear in the lattice theory. The energy eigenstate \Bq{x)) can be thought of as composed 
of a static quark fixed to the position x together with a light-quark bound state centered 
at X. The tilde notation is chosen to suggest the relation of these new states to our earlier 
energy and momentum eigenstates, \B{'p)v), which are standard superpositions of these 
Nh{x) eigenstates: 

m'p)y) = Y,j^\BM)- (CIO) 



Using these Nh{x) eigenstates and the conservation of Nh.{x') for all x' G ^ it is now 
straightforward to derive Eq. (!72l) . Substituting a sum over the complete set of degen- 
erate ground states \Bq{x)) into the correlation function C^^{t,to) which appears in the 
denominator of Eq. (172]) and is defined in Eq. (1671) we find: 

C^^(t,to) = E ((0|Mx)7o75(E ^(y))\B,m (Cll) 

xeAV I yeAV 

•(5,(x)|( J2 ^(r))7o75/^(x)|0)|e— 

y'£AV J 

where we have used the conservation of local heavy-quark number to keep only terms where 
all of the heavy quarks are located at the same spatial position x. 
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Similarly, we can evaluate Co-(tf,t,to) which appears in the numerator of Eq. (172|) and 
is defined in Eq. fl69|) : 



Cg^(t,,t,to) = V ■! {0|( V mt)h,,7M!.t)\B,m (C12) 



xeAV I 



y&AV 



■{B,{x)\0,{x)\B,{x)) ■ ^ qiy',t))^ol,hix,t)\0) U-'^-^*/"*") 

= C^^{tf,to)(B,i6)\0,{0)\B,{0)), (C13) 

where translational invariance of the matrix element {Bq{x)\Oi{x)\Bq{x)) has been used to 
remove it from the sum over x causing that sum to assume the same form which appears 
in the box-box correlator C^^{t,to) given in Eq. ( ICllI) . after an application of charge con- 
jugation symmetry. Equation ( 1721) is then easily recognized from ratio of Eqs. ( ]C12|) and 

dcn]): 

Mo, = ^(Bq{p = 0)\OS)\Bqip = 0)) (C14) 

= 2L'(Bq{p = 0)vm6)\Bq{p = 0)v) (C15) 

= 2(^,(0)|O,(0)|S,(0)) (C16) 

= 9. ^' ' C18 

where Eq. fICSp . relating our two normalization conventions for momentum eigenstates, has 
been used to obtain the second equation. 

This use of localized sources while making the static approximation may become more 
important as larger spatial volumes are used and the overlap between the translationally 
invariant wall sources and the physical states of interest becomes smaller. Of course, this 
same method can be used for localized sources with different spatial distributions such as 
Gaussian or atomic wave functions. 



Appendix D: Chiral perturbation theory for decay constants and mixing parameters 

In this section we present the NLO HMxPT expressions needed to extrapolate A'^^ = 2 + 1 
domain-wall lattice data for heavy-light meson decay constants and mixing parameters to 
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the physical quark masses and the continuum. Although we label the formulae with the 
subscript "5" , we note that these functions can also be used to extrapolate _D-meson decay 
constants and mixing matrix elements with the caveat that the low-energy constants (except 
for the light-light meson tree-level parameters / and B) are different for the case of 5-mesons 
and D-mesons. We first show the SU (3) HMxPT formulae in Sec. ID 1( we then take the 
appropriate limits of the SU (3) expressions to obtain those in SU{2) HMxPT in Sec. ID 2[ 



1. SU{3) HMxPT expressions 



The tree-level mass-squared of a meson composed of two domain- wall valence quarks with 
flavors X and y is 



m^y = Bim^ + my + 2m^ 



(Dl) 



where B is a. continuum low-energy constant and m^cs is the residual quark mass. 

The NLO result for = fn^-J mR^ in the partially-quenched domain- wall theory with 
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2-1-1 flavors of sea quarks is 

1 1 + ^g^B. 



82|: 



1 



1 1 + ^9%*B. 



f=l,l,h 



6 



je{A/x} 



^(i?fk{Mx};{/i}))^K) 



(D2) 



where / ~ 130.4 MeV is the tree-level pion decay constant. The NLO expression for Mb^ 



8/3 TUB^fB Bb^ is similar 



M, 



1 + ^9l*Bn 



2 \ 1 '^dB'B-K 



-f 1 1 u J 



d 



2A7i^P 



je{Aix} ^ 



+ dsea(2m/ + rrih) + d^aim^ + daO? 



(D3) 



In both the decay constant and the mixing matrix element, the only effect of the nonzero 
lattice spacing is a new analytic term proportional to a?. These results agree with the 
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continuum calculation of Sharpe and Zhang in the limit of three degenerate sea quarks 23 |. 
We note that one is free to multiply the above expressions by arbitrary powers of the heavy- 
light meson mass m^^ without modifying the chiral logarithms. This is because the chiral 
logarithm contributions to heavy-light meson masses are suppressed by 1 / m^, and therefore 
of higher-order than we consider jSJ] . 

In the above expressions for and Mb^, the functions i and i are one-loop chiral 
logarithms: 

,2 



m In — 



1 . 



(D4) 
(D5) 



The function is due to single poles in the partially-quenched propagator: 



and the sets of flavor- singlet masses that appear in the residue functions are 

{Mx} = {m\,m\} , 
where is the mass-squared of a meson composed of two x valence quarks, 



(D6) 



L{H) 



(D7) 
(D8) 

is the 



mass-squared of a meson composed of two lih) sea quarks, and = (m| -|- 2m^)/3 for 
2-1-1 flavors of sea quarks. 

For completeness, we also include the expressions for the decay constant and mixing ma- 
trix element at the unitary points. In the full QCD and isospin limits, the above expression 
for $a; becomes 



1 



167r2/^ 



1 + "igl^B. 



(D9) 



1 - 



167r2/2 



1 + 'igl*B. 



+ Csca(2m; + nih) + C^i^mh + CaO^ 



(DIO) 
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where = m| and mj^ = {m\ + m?^) /2. Similarly, the expression for Mx becomes 

^ [3^(m^) - ^(mj)] + 4ea(2mz + m,,) + d^^mi + 4a' > , (Dll) 



487r2/2 



Mb, = Poll- ,rJT MO + ^K)] - 



167r2/2 L V ii/ ' V i/J iQ^2p 



+ ^^^2^-2 - 2^(m2)] + 4ca(2mi + nih) + d^^imt + daO? ) . (D12) 



The expressions for and agree with those derived by Goity in Ref. [J 

One can account for lattice finite volume effects at NLO in xPT by turnin g th e one-loop 
integrals to sums. This yields an additive correction to the chiral logarithms [86|: 

«K)=m^(l„^ + C-(™L)). ir(™L) = -l-^M^, (D13) 



- L^ + l] +53^^(mL), 53^^(mL) =25^i^o(|r1mL), (D14) 

V X / i^^Q 

where 5f^{mL) is the finite volume correction to the infinite volume result and Kq and Ki 
are modified Bessel functions of imaginary argument. 



2. SU{2) HMxPT expressions 

The NLO SU (2) HMxPT expressions can easily be obtained from the SU (3) results in 
the previous subsection by integrating out the strange valence and sea quarks. After this 
procedure, however, the expressions for the decay constant and mixing matrix element differ 
for i?(i-type mesons and -B^-type mesons. 

First we consider i?(i-type mesons composed of a 6-quark and a light valence quark with 
mass TJix- In this case, we take the limits of Eqs. f lD2p and (IDSp assuming 

(D15) 

vfih rrih 

The resulting expression for at NLO in the partially-quenched domain-wall theory with 
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two degenerate light sea quarks is: 




(47r/(2)^ 



(4vr/(2))2 



{ml - m\Y{m\) + i{m\) 



(D16) 



where we use the superscript "(2)" to distinguish the SU{2) low-energy constants from their 
SU (3) analogs in the previous section. The NLO expression for Mb^ is similar: 

l-3(9g.V)^ 



M. 



(2) 



(4vr/(2))^ 



2£(m: 



(4vr/(2))^ 



+ 



(mi-m^K(m^)+£(m^; 



+ d^^jni + df},m, + d^^^a^\ . (D17) 



(4vr/(2))^ 

These results agree with the continuum partially-quenched calculation of Sharpe and Zhang 

The unitary QCD expressions can easily be obtained by the 



in the limit a — )■ 



replacement m^; mi. In the SU{2) theory, the effects of the dynamical strange quark are 
fully contained in the values of the low-energy constants, e.g. 0q {m^hj- For simplicity of 
notation, however, we suppress the functional dependence of the coefficients on m^. 

Next we consider iJ^-type mesons composed of a 6-quark and a heavy valence quark 
with mass my. Because the SU{2) chiral effective Lagrangian includes only two light quark 
flavors, this requires an extension of SU (2) xPT to the kaon sector, and the resulting theory 
is sometimes called Kaon SU (2) Chiral Perturbation Theory (K^PT) 
we take the limits of Eqs. ( ID2p and flD3l) assuming 
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87|. In this case. 



mi mi 

— ,— <!. 

my mt 

The resulting NLO expressions for $5 and Mb are: 



(D18) 



By 



l + ciimi + c):'a 



Wn2 



(D19) 



By 



Pi'Hl + diimi + S:^a' 



(D20) 



where we use the superscript "(s)" to distinguish the coefficients from those in Eqs. (1D16P 
and (ID17p . Because the valence quark has been integrated out, there are no longer any 
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chiral logarithms. The effects of both the valence and sea strange quarks are encapsulated 
in the values of the low-energy constants, e.g. (j)Q\my, nih). 
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